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We have studied the effect of impurities on the criticai temperature for a model of YBCO involving 
pairing both in the Cu02 planes and in the CuO chains. In this model pairing in the planes is due 
to phonons, while Coulomb repulsion induces in the chains an order parameter with opposite sign. 
Due to the anticrossing produced by hybridization between planes and chains, the order parameter 
changes sign on a single sheet of the Fermi surface resulting in nodes in the gap. We find that, 
in our model, the critical temperature is much less sensitive to impurities than in standard d-wave 
models. One reason is that impurities produce essentially plane-plane and chain-chain scattering, 
which does not affect the critical temperature. T c is reduced by the scattering between parts of 
the Fermi surface which have opposite signs for the order parameter, just as in standard d-wave. 
In our model this is due to plane-chain scattering. We have found that this scattering, whatever 
its origin, will be smaller by a factor of order t /Ef ( that is hybridization coupling over Fermi 
energy ) compared to plane-plane and chain-chain scattering. Accordingly the sensitivity of T c to 
impurities in our model is reduced by a similar factor, compared to the d-wave situation. In the 
specific case which we have studied in details and which reduces to the two-band model, we have 
found a further reduction of the sensitivity of T c to impurities with a behaviour which can vary 
continuously from s-wave like to d-wave like depending on the parameters. We expect a similar 
behaviour and reduction to occur in the general case. 

PACS numbers : 74.20.Fg, 74.72.Bk, 74.25. Jb 



I. INTRODUCTION 

The mechanism responsible for high T c superconductivity is still highly controversial jjj . While a large part of the 
theoretical effort is based on the hypothesis that Cooper pairs are formed in high T c compounds, there is no agreement 
on the physical origin of the pairing interaction nor on the symmetry of the pair wavefunction. These two questions are 
actually intimately related. Indeed repulsive interactions, such as produced by spin fluctuations Q|, require an order 
parameter which changes sign on the Fermi surface in order to produce pairs. This leads, in the simplest hypothesis, 
to pairing with d-wave symmetry if we assume singlet pairing. On the other hand purely attractive interactions lead 
to pairing with s-wave symmetry since any change of sign is unfavorable in this case. Therefore an experimental 
determination of the order parameter symmetry should greatly help to identify the physical interaction responsible 
for pair formation: although it would not be enough to provide a unique identification, it would strongly narrow the 
remaining possibilities. For this reason a large part of the recent experimental work has been aimed toward providing 
a clear signature for the symmetry of the order parameter. 

Quite surprisingly recent experiments on Y BCL2CU3O7 designed with this purpose of identifying the symmetry have 
given clear cut, but contradictory answers [Q. Indeed the observation of a sizeable Josephson current |}| in a c-axis 
tunnelling junction between Y Ba^Cu^O-; and Pb is quite difficult to reconcile with a pure d-wave symmetry while 
it is in full agreement with s-wave symmetry Similarly the fact that there is no angular dependence in the critical 
current of YBCO - YBCO grain boundary junctions in the a-b plane 0] goes clearly in the direction of an s-wave 
interpretation. On the other hand a number of experiments are in favor of a d-wave symmetry. Many experiments, 
including tunnelling, NMR, Raman scattering, photoemission, penetration depth have shown the existence of low 
energy excited states. Actually these experiments are compatible with a strongly anisotropic s-wave order parameter. 
Or more simply one may look for extrinsic effects and wonder if these states do not arise from surface effects or defects, 
which would provide quasi- normal regions. However the existence of a linear T dependence of the penetration depth 
over a large range of temperature, simultaneously in cristals and films, with the same slope makes an extrinsic 
interpretation for all these experiments unlikely ( while this kind of explanation may very well be valid for some 
of them ). Moreover some experiments specifically designed to check if the order parameter changes sign over the 
Fermi surface have given positive answers. These are the corner SQUID experiments Q which give a clear indication 
for a change of sign of the order parameter between the a and the b axis, and the observation of a spontaneous 
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magnetization corresponding to a half magnetic flux quantum in 3 grain-boundary Josephson junctions m which 
implies a n shift, in clear agreement with d-wave symmetry. 

On the other hand the d-wave interpretation is not free of problems. For example recent experiments show that 
some thermodynamical superconducting properties are markedly anisotropic ( the a and b-axis results are different 
). Indeed the penetration depth in good YBCO cristals display a strong anisotropy of the penetration depth ||, the 
specific heat anomaly at T c of the parent compound LuBa2Cu^Oi ( with same T c as YBCO) has a marked anisotropy 
as a function of the orientation of an applied magnetic field ft is difficult to ascribe these anisotropies to the 
weak orthorhombic distorsion, and it is more likely that the CuO chains play a significant role in the superconducting 
properties. One of the most conspicuous problem of the d-wave interpretation is the weak sensitivity of the critical 
temperature of Y Ba2Gu^Oi to the presence of impurities. Indeed any kind of impurities, whether magnetic or not, 
produces in d-wave superconductors JlO|— p~2[ an effect analogous to pair-breaking by magnetic impurities in standard 
s-wave superconductors |l3| . In particular the critical temperature pO| , [Ll| decreases rapidly with increasing impurity 
concentration following the Abrikosov - Gorkov law [ fl3[ , and superconductivity disappears at a critical concentration. 
In contrast all samples of Y BaiCu-sOi seem to have a T c around 90 K. It is difficult to believe that all samples ( 
including the earlier ones ) are clean enough to affect only weakly the critical temperature, whereas all microscopic 
studies show that there are always more structural defects than what is generally admitted. Moreover ion (l4|] or 
electron fl5f| irradiation experiments have shown a rather weak sensitivity of the critical temperature of Y BCL2CU3O7 
on the inverse lifetime deduced from resistivity measurements [OJ. Actually Zn impurities are known [ p^|Jl7[ to have 
a depressing effect on T c , but this can be interpreted as a standard pair-breaking effect since the environment of a Zn 
impurity, located in the CuOi planes, is known to be magnetic ]lS| ]. 

In order to solve these contradictions we have proposed recently JljJ for Y Ba 2 Cu^Oj a model which mixes s-wave 
and d-wave features. In our model , in addition to the CuOi planes, the CuO chains play an essential role. The 
pairing interaction within the planes is attractive ( it can be for example produced by phonons ). On the other hand 
the pairing interaction between planes and chains is repulsive ( it can be produced by Coulomb interaction ). In this 
way the order parameter has opposite signs on the planes and on the chains. Moreover we include the hybridization 
between planes and chains, which corresponds physically to take into account the possibility for an electron to jump 
from planes to chains or vice- versa. Naturally the coupling responsible for this hybridization is fairly small, but it is a 
well known feature of all band structure calculations [|20f| . It is of importance only when the plane and the chain band 
intersect. In this case it leads to an anticrossing in the dispersion relations, and similarly to an anticrossing of the Fermi 
surfaces, wherever the (uncoupled) pieces of the Fermi surface related to plane and chain cross. As a result, when we 
move on a given sheet of the Fermi surface, we go from a part which corresponds physically to a plane electron, to 
a part which corresponds physically to a chain electron. Since the order parameter has opposite signs for plane and 
chain electrons, this implies that the order parameter changes sign on a given sheet of the Fermi surface and therefore 
has nodes on this sheet by continuity. Therefore our model provides an order parameter which is quite analogous to 
a d-wave order parameter and it can in this way explain [fl9[ all the experiments in favor of d-wave symmetry. On the 
other hand it does not have a d-wave symmetry since the nodes, which occur at the intersection between plane and 
chain bands, have no reason to satisfy k x = k y , and the average of the order parameter has no symmetry reason to 
be zero. Hence there is for example no problem with the nonzero Josephson current in Y BCL2CU3O7 - Pb junctions 
along the c - axis || . Moreover since the attractive in-plane interaction and the repulsive interaction between plane 
and chain help each other, there is no problem to explain the high value of the critical temperature of Y Ba 2 Cu^Oj 
plfl . Naturally our model is specific of YBC12CU3O7 , but we can think that it can be generalized to other compounds 
where the role of the chains can be played by other parts of the structure, such as the BiO planes in BSSCO. On the 
other hand there is no possibility of this kind in LSCO and accordingly we do not expect experiments to display in 
this compound the same physical features as in YBCO. 



Our model has common features with many other models. As a two band model it is quite similar to the two 
band models introduced by Suhl, Matthias and Walker and others fj2|-p6f to describe superconductivity in transition 
metals. The possibility of an interband repulsive interaction in two band models was already introduced by Kondo 
pi| . More recently a two band S-N model has been introduced by Abrikosov and Klemm ^Z7j to account for Raman 
scattering data. The idea of an interband repulsion has been put forward by various groups p8|-^l| recently in the 
context of high T c superconductivity in order to show that experiments displaying a change of sign of the order 
parameter did not necessarily imply a spin fluctuation mechanism. Our model is also similar to the one proposed by 
Abrikosov |32|, where there are attractive and repulsive interactions within a single band, leading to a change of sign 
of the order parameter within this band. Finally, as noted above, chain-plane hybridization is a standard feature of 
band structure calculations pG] and there have been various suggestions of the importance of the chains in the physics 
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of Y ' Ba-iCuzO-i , whether for intrinsic or for extrinsic reasons |33|-|35[|, these models allowing for the electron to jump 
between planes and chains. 

In this paper we consider the effect of impurities on the critical temperature in our model. We will actually restrict 
ourselves to non magnetic impurities. Indeed magnetic impurities are easily included, but they will naturally lead 
to pair-breaking and produce a T c following the AG law as in other models. Hence magnetic impurities can not 
discriminate between various models and we ignore them for simplicity. The conclusion of our study is that, in our 
model, for generic parameters, the critical temperature is much less sensitive to impurities than in standard d-wave 
models. Therefore we might say that, with respect to the sensitivity of the critical temperature to non magnetic 
impurities, our model behaves as a "weak" d-wave model. This happens for two independent reasons. First the 
reduction of T c is due to plane-chain scattering, which is weak compared to plane-plane scattering. Next the fact 
that we have a two band model provides further possibilities for a weak impurity sensitivity. In section II we present 
our model and we calculate the critical temperature of the clean superconductor as a first step toward the calculation 
of the impurity dependence, which is dealt with in section III. Finally section IV is devoted to a comparison with 
experimental results and to our conclusion. 



II. THE MODEL 



Let us first define specifically our model. We will describe the motion of an electron in planes and chains by the 
following Hamiltonian: 

H = ^ £kc~lc k + ^ £ 'kdtdk + ^ t k c t d k + h.c. (1) 

where ct and are creation operators in the plane and in the chain band respectively. The first term corresponds 
to an isolated plane with dispersion relation e^, the second one to isolated chains with dispersion relation e' k and the 
last term describes hopping between planes and chains. Actually this independent electron Hamiltonian does not 
correspond precisely to the situation found in Y BazCu^Oi . Indeed YBaiCuj l O'j is build by stacking up sets made 
of two CuOi planes and one CuO chains plane. Therefore a better description is obtained by the Hamiltonian: 

H = J2 ec l.n c l,n + E £c 2,n C 2,n + E tp( C l,n C 2,n + h.C.) 
n n n 

+ E e' d n d n + E tcict^dn + h.c.) + E t c (c^ n d n+1 + h.c.) 

n n n 

where all quantities are understood to depend on k x and k y and summations run also over k XtV . The indices 1 and 2 
number the CuOi planes and the index n numbers the stacks. Introducing the even and odd plane band operators 
c± by ci.2 = ( c + ± c_ )/\/2 , and taking the Fourier transform in the z direction, we obtain: 

#o = E( e + tp) c +,k c +,k + E( e - t P )ct )k C-,k+ 

( 3 ) 

+ Y,e'd+d k + E*+(c|,fc d fc + h.c.) + + h.c.) 

where a dependence and summation on k z is now understood, and t + = t c y/2 cos(fc z c/2) and t _ = t c \j1 sin(fc z c/2) 
, with c being the size of the unit cell along the z direction. Band structure calculations J2(| give a crossing of the 
chain band Fermi surface with the odd plane band Fermi surface, whereas there is no crossing with the even plane 
band. Therefore the even plane band will not play an interesting physical role in our model and we will forget it 
for simplicity, but it should clearly be retained in a very realistic calculation. This leads us back to the Hamiltonian 
Eq.(l) where the plane band is actually the odd plane band. The hopping term has a k z dependence as given by t _. 
However, if the even plane band happens to cross the chain band, we would have no problem in extending our model. 
This would be particularly easy in the two extreme cases where the even and odd plane bands are nearly degenerate 
(t p <<t c )or well separated ( t p >> t c ) since we would be led again to the Hamiltonian Eq.(l). 

With respect to pairing interactions, as we have already indicated, we take an attractive pairing interaction in the 
plane which can be for example due to phonon exchange. On the other hand we assume a repulsive pairing interaction 
between plane and chain. A natural physical origin for this term is the Coulomb repulsion, since screening is certainly 
not very efficient in YBCO because it is a weak metal. Moreover screening is probably very ineffective along the z 
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direction since hopping in this direction is small. We note however that hopping is physically necessary in order to 
have a pairing term between planes and chains. Finally there is little experimental evidence for an attractive pairing 
interaction in the chains. However we will assume its existence for the sake of generality. We include also the standard 
intraband Coulomb repulsion. 

Let us first consider the critical temperature of the clean superconductor. In a first step we ignore the hopping 
term t The Eliashbcrg equations at the critical temperature are : 

A n Z n = ^VVmpi - ttT/x V ^ - nTk V ^ (4) 

A' n Z' n =nTj2 A'n-m - V ^ _ ^ Tfc ' g ^ (5) 

Here A„ and Z n are the order parameter and the renormalization function at the Matsubara frequency u> n . Primed 
quantities refer to the chains while unprimed ones correspond to the planes; A„ and A^ are the effective frequency 
dependent interactions due to phonon exchange while /i, // and k, k' are the intraband and interband Coulomb 
repulsion with cut-off uj c of order of the Fermi energy. While the interactions that we consider here have certainly 
a wavevector dependence in the real superconductor, we have no precise idea of what they are nor do we have any 
physical reason to believe that this dependence is strong. Therefore, for the sake of simplicity and clarity, we have 
taken isotropic interactions within planes and chains. We can now proceed to the usual reduction of the cut-off from 
lj c to a frequency lojj a few times a typical phonon frequency. For lu„ > cup , the order parameters A n and A' n become 
constants Aqo and AJ^ while Z n , Z' n w 1. These constants are obtained from Eq.(4) and (5) as : 

(l + M r)A 00 +A;rA / 00 =-7rT/iV^-7rTfcVf^ (6) 



(i + A)a'oo + fcVAoc = -tttv Y,rn~ * Tk ' E ( 7 ) 



where r = In (to c / lod)- When these results are carried into Eq.(4) and (5), one obtains the same equations, except 
that the cut-off is now u>£> and /i, // and k, k' are replaced by renormalized Coulomb interactions fi*, fi'* and k*, k'* 
given by : 

* = /i + (/i// -kk')r = k 

^ (1 + jur)(l + fi'r) - kk'r 2 (1 + /xr)(l + n'r) - kk'r 2 1 ' 

and similar expressions for //* and k'* obtained by exchanging primed and unprimed quantities. Naturally one 
recovers the standard Coulomb pseudopotential when there is no interband interaction. It seems safe to assume 
- kk' > , otherwise we would have pairing from 's-wave' Coulomb interaction alone (when pairing from Coulomb 
interaction is considered, this is usely for higher 'angular momenta' ). In this case renormalization does not change the 
sign of the various quantities. We see that, in addition to its direct effect, interband interaction decreases the effective 
intraband repulsion, leading to an increase of the critical temperature as it should fl36| . In the following we will 
restrict ourselves to the weak coupling limit where A„ and A' n can be taken as constants and mass renormalization is 
negligible. Indeed, as it will be clear in the following, we do not expect strong coupling effects to modify qualitatively 
our conclusions. In this case T c depends only on the combinations A* = Ao - /i* and A'* = Aq - (J, 1 *, in addition to k* 
and k'*. From now on we will omit the star for all the renormalized coupling constants. Then, if we let x = In ( fi / 
T c ) where fi is a typical phonon frequency ( or rather 1/4 of a typical phonon frequency j3?J ) , T c is obtained from : 

(aT 1 - AXaT 1 - A') = kk' (9) 

which is the standard two band result |22]. Although one expects physically that Coulomb interaction reduces the 
isotope effect, this is not obvious from Eq.(8) and (9). This has been shown to be true by Kondo p4| in some limiting 
cases. We show in the Appendix that this property holds quite generally. 

As a preliminary to the calculation with impurities, let us finally consider how the above calculation of the critical 
temperature is modified when we take into account the hybridization t k between planes and chain. As indicated 
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above, we take isotropic interactions and order parameters in the planes and in the chains. Corresponding to the 
effective interaction Hamiltonian : 

H int = -g^2 c i> c ~t k ,c- k ck + K^d+dt k ,c- k c k + h.c. - g'^d+dt k ,d-kd k (10) 
the order parameter satisfy : 

A = -gJ2< c -kC k >+Kj2< d -kdk> (11) 
k k 

A' = K < c -kc k > -g' < d - kdk > ( 12 ) 
k k 

However because of the hybridization, the plane and chain operators Ck and dk do not correspond to the eigenstates 
of H . The unitary transformation : 

c% = (Hjlj (13) 

diagonalizes the Hamiltonian into H — S ej 7+ 7$ where we have set for convenience Ci = c and c 2 = d . The energies 
ei and e 2 of the hybridized bands are given by 2 ei i2 = e + e' ± [ ( e - e' ) 2 + 4 t 2 J 1 / 2 and an = a 2 2 = [(ei - e')/( e i 
- e 2)] 1 ^ 2 = cos 6* , a 2 i = - ai2 = [(ei - e)/(ei - e 2 )] 1 / 2 = sin 6*, with again the k dependence understood everywhere. 
When we carry out this transformation in the mean field Hamiltonian : 

H = H + A £ c+c+ fc + A' <dt k + h.c. (14) 

k k 

we will obtain interband pairing terms, such as 7_fc l i jk.2 , coupling the hybridized bands. However we will consider 
that t k is large enough so that these bands are well separated. Specifically this means that we assume ujd « t k ( 
otherwise we should make a much more careful strong coupling treatment, which would probably not bring anything 
new qualitatively). In this case the above pairing terms will be negligible because we can not have two electrons (k,l) 
and (-k,2) at the Fermi surface with opposite wavevectors, but belonging to different bands ( we will justify more 
specifically this approximation in the next section ). Therefore the transformation Eq.(13) gives : 

H = H Q +Y, Kkltklt-k + J2 5 ^kltkll-k + h - c - (15) 

k k 

with 5\ t k — A cos 2 k + A' sin 2 k and 5 2 ,fe = A sin 2 6 k + A' cos 2 6>fc. Naturally this band diagonal expression for 
the Hamiltonian is obvious once the interband pairing terms are neglected, and we could have written it immediately. 
Our essential point is that, starting from isotropic interactions in planes and chains, we obtain a specific anisotropy 
for the order parameter 5i tk and <5 2i fc . In particular we obtain nodes at the Fermi surface since we have managed to 
have A and A' with opposite signs, and cos 2 k goes from to 1 when we move at the Fermi surface of a given band. 
Therefore we have an order parameter which is d-wave like, although we have assumed an attractive pairing in the 
planes. 

We have then as usual < 7 7 j > = - (<5j/2ej) tanh(/?ej/2). From Eq.(ll),(12) and (13) this leads to: 

A = V(gcos 2 fl fe - Ksin 2 9 k )^- tanh(^) + (gsm 2 e k - Kcos 2 e k )p±- tanh(^£) (16) 

2ei >k 2 2e 2 .k 2 

A' = V(.9'sin 2 fl fe - Kcos 2 k )^- tanh(^M) + (g'cos 2 9 k Ksin 2 6 k )^- tanh(^) (17) 

^ tel,k & ^2,k ^ 

Since the bands are well separated, the k summations will be around the Fermi surface of each band with cut-off u>u 
for e^fc . The integration perpendicular to the Fermi surface will give the standard factor x = In ( 1.13 (3 u>d) and we 
are left with summations along the Fermi surface: 
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A = zy^(gcos 2 6> fc - Ksm 2 9 k )5 lyk 8(e 1 , k ) + (gsm 2 9 k - Kcos 2 9 k )5 2tk 5(e 2tk ) 

k 



(18) 



A' = 



a; 



^(g'sin 2 (9 fe - Kcos 2 9 k )6 hk S(e 1 . k ) + (g'cos 2 9 k - Ksm 2 9 k )8 2 . k S(e 2 . k ) 



(19) 



k 



where the order parameter Si <k should not be confused with the Dirac function 8(e). For e% <k = we have cos 2 9 = 
e'/(e + e') and cos 2 9 = e/(e + e') for e 2 , k = 0. Moreover for most of the Fermi surface we have e << e' or e' << 
e leading to cos 2 9 = or 1. This corresponds merely to the result one obtains without hybridization. Subtracting 
from Eq.(18) and (19) their counterpart without hybridization, we are left with summations over quantities which are 
essentially nonzero only in the anticrossing region. In this small region it is convenient to take, at fixed k z , e and e as 
new variables instead of k x and k y . The Jacobian of the transformation is J = d (e, e')/d (k x , k y ) = | v p x v c | where 
v p and v c are the Fermi velocities of the plane and chain band at the crossing point. We are then led to evaluate : 



with /(e,e') = e 2 or e' 2 or — ee'. This is easily done by taking e/t k and e'/t k as new variables and extending the 
boundaries of the integrals to infinity. The result is the same for the three integrals, namely - it \t k |. When this 
result is carried into Eq.(18) and (19) we find naturally the first order correction to the coupling constants A = N p g 
, A' = N c g' , k = N c K and k' = N p K , where N p and N c are the total density of states of the (uncoupled) plane 
and chain bands. The hybridization changes these coupling constants respectively into A - N t (g+K), A' - N t (g'+K) 
, k - Nt (g+K) and k' - Nt (g'+K). Here Nt is an effective density of states due to the hybridization and given by 
N — i /(irJc) where i is the average of \t k | over k z (we have taken into account that there are 4 crossings in the 
Brillouin zone ). Naturally since we have Nt << N p , N c in our model these changes are small. 



We consider now the effect on the critical temperature of impurities located in the planes and in the chains. We 
will assume that the impurity potential scatters electrons within the planes or within the chains, but that it does not 
scatter them from plane to chain. Our physical motivation is naturally that planes and chains are physically well 
separated, which corresponds to the fact that the hopping term t k is small. Therefore the potential for scattering 
from plane to chain is reduced by a factor of order t k / to compared to the intraplane or intrachain scattering ( 
where to is a typical hopping term within planes or chains ) . This physical point is essential for the validity of our 
result. Including this small plane-chain scattering does not make any problem, but for sake of clarity and simplicity 
we will not do it explicitly and we will rather come back to this point after we have obtained the effect of intraplane 
and intrachain scattering. We will also take an isotropic scattering within plane and chain, consistently with what 
we have done for the interactions. Here again we will come back later on to the effect of a possible anisotropy. This 
leads us to the following impurity potential : 



We treat this potential within the Born approximation since, for the calculation of the critical temperature, making 
use of the T-matrix approximation is merely equivalent to renormalize the scattering potential. There is however 
another effect produced by going beyond the Born approximation, which will be considered at the end of this section. 

Within the Born approximation and after impurity averaging, only contributions to the self-energy corresponding 
to scattering twice on the same impurity are retained p8[ . Since the potential Eq.(20) does not scatter electrons from 
plane to chain, the impurities contribution to the self-energy has only components within planes or chains. They are 
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(20) 



k,k> 



k,k' 



given by [|9|: 




(21) 



k 



G 



where n\ and n 2 are the number of impurities per unit volume in the planes and in the chains, Tj are the Pauli matrices 
in Nambu space and G is the temperature Green's function in plane-chain representation. The Green's function is 
related to the self-energy by : 

G" 1 = Gq 1 - X (22) 
where Go is the Green's function in the absence of the impurity self-energy . Explicitly we set : 

G"- 1 = (iCli - SiT 3 - AiTi)5ij - tT 3 a x , io (23) 



G o,lj = ( iuj ~ e » r 3 _ - tT 3 a Xtij (24) 

where u> = (2n+ 1) it T is the Matsubara frequency ( we omit systematically for clarity the index n in the Matsubara 
frequency and all frequency dependent quantities ), a x is the Pauli matrix in plane-chain space and Aj is the off- 
diagonal self-energy due to the pairing interaction ( we have set ei = e and e 2 = e' )• In order to invert G _1 it is 
convenient to go the hybridized band representation, or more precisely to make the transformation which diagonalizes 
G _1 when the off-diagonal part Aj of the self-energy is zero. For Gg -1 the transformation matrix is just To where 
dij is the transformation considered in the preceding section. In the present case we have merely to formally generalize 
it by replacing ti by - i Oj for the particule-particule part in Nambu space and by ej + i fij in the hole- hole part. We 
keep the same notation for this transformation, but the Ojj are now complex. The complete transformation matrix is 
now Aij — ( dij + a*j ) tq /2 + ( - a*^ ) T3 /2 . Actually we have X , A << t since we assume that the hybridized 
bands are well separated ( X or A are typically of the same order of magnitude in the important cases as it will be 
clear below ). Moreover as usual we will be only concerned with Matsubara frequencies u at most of order S and A 
and therefore negligible compared to t . This implies that, when it multiplies a quantity of order X or A , we can 
take for its zero order value, namely the real value found in the preceding section. This leads to the following 
expression g^ 1 = l A G _1 A for the full Green's function in the hybridized representation ( *A is the transpose of A ) 

glj 1 = - e*i) T o + (e» + e* )r 3 ]% - a ki A k akjTi (25) 

where : 

2ei, 2 = ei + e 2 - + f2 2 ) ± [(£1 - £2 - i&i + i&if + 4t 2 ] 1/2 w 2e h2 - 2ia>i, 2 (26) 
and we have set : 

2u, li2 = n 1 + n 2 ± ("i-"')( g i- e ') (27) 

Now the dominant contributions to the summations of G over the wavevector k will come from the vicinity of the 
Fermi surface of the hybridized bands. In this case either e\ is large and all the other matrix elements of g^ 1 are 
small ( of order A typically ), or e 2 is large and the other elements are small. It is then easy to see that when e 2 is 
large, we obtain g by simply inverting the (1,1) block of g -1 , that is have gij w [(.9 ^n] 1 <5j,i- The other terms 
of g^ are smaller by a factor of order A over energy separation of the hybridized bands ( for fixed k ) , that is a factor 
of order A/t or less . We have a similar result when e\ is large. Taken together this means that we can take g^ as 
block diagonal, that is — gn Sij with : 

e l e* + Sf V Oi -ei / 

where Si — af^ Aj , that is explicitly 61 — Ai cos 2 6 + A 2 sin 2 6 and <S 2 = Ai sin 2 8 + A 2 cos 2 6. 

We can now calculate the impurity self-energy S . Since g^ is block diagonal we have Gu — a 2 gjj . When we 
carry this result into Eq.(21), we can perform the integration perpendicularly to the Fermi surface (that is integrate 
over e\ or e 2 ) since it converges within an energy range of order A. We are then left with summations along the 
Fermi surface, just as in the calculation of the critical temperature in the preceding section. This leads to : 
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k 



(29) 



where we have taken into account that Si << uii since we are only interested in the calculation of the critical 
temperature. We can then proceed as in the preceding section : for most of the Fermi surface we have e\ « e 2 or 
e 2 << ei leading to afj = or 1. This corresponds to the result without hybridization which is : 



-Ti(iT signuj -n) 



(30) 



where I\ = tt rn AT; Uf with N\ = N p and N 2 = N c . When we calculate Si - S° we are left with summations over 
nonzero quantities only in the anticrossing region. We take then (ei - e 2 ) / t and (ei + e 2 )/ 2t as new variables and 
extend the boundaries of the integrals to infinity. The resulting expression simplifies under exchanging ei and e 2 • 
We obtain that the diagonal parts of Si and S° are equal which implies Hi = oj + Ti sign uj . On the other hand we 
find from the off-diagonal part : 



Ai 



^ Ai N t 2 f 1 du 

Al + ri — ^ + r 1 -^- / — 

Oi ^iWo (l_ u 2)3/2 



(1 + w,) 2 Ax + (1 - u 2 )A 2 



(i + «)fii + (i-«)n 2 



+ {u — > -ti} 



2Ai, 



Oi 



(31) 



where u = (ei - e 2 )/[(ei - e 2 ) 2 + 4 t 2 ] 1/<2 , and a similar result for A 2 



This equation displays explicitly an essential point for our result : the correction to the self-energy due to hy- 
bridization is small, of order N t / Ni t2 , compared to the self-energy in the absence of hybridization. This correction 
corresponds physically to the scattering from plane to chain induced by hybridization. However, without hybridiza- 
tion, we are in the standard s-wave situation and impurities do not change the critical temperature. The only effect 
on T c will come from the correction due to hybridization. Since this correction is small, we expect the change of T c 
to be small with respect to what one would get with a simple d-wave order parameter. Although the integration 
in Eq.(31) can be performed analytically in the general case, the result is not simple. Therefore, in order to obtain 
reasonably simple calculations, we will continue our quantitative investigation only in the particular case where Ti = 
r 2 = r, which implies Cli = Q 2 = Q . Actually this equality is always valid when the Matsubara frequency to is large 
enough compared to Ti and T 2 . This is in particular always valid when the impurity concentration is small enough 
so that i\ 2 << T c . Therefore the case we consider is quite reasonable. 

In this case the equations for Ai and A 2 simplify into : 



Ai = Ai + r 



Ai 



N t A 2 - Ax 



AT, 



(32) 



a 2 = a 2 + r 



A, 



,N t A\ — A 2 



A^ 2 



(33) 



We have finally to calculate from Eq.(12) and (13) the self-energies Ai and A 2 due to the pairing interactions. We 
have already obtained the off-diagonal part of the Green's function, and the calculation is essentially the same as for 
the off-diagonal part of the impurity self-energy. One finds : 



Ai = A^tV[^ + 

o 

77. ^ M 



N t A 2 — Ai AT jz ^^ r A 2 N t A1-A2 
M 1 ~ NiK-kT ZJ.tzt + tt : ] 



(34) 



a 2 = W£ A + - w£[^ + 



(35) 
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However these equations merely mean that we have to calculate the critical temperature with the coupling constants 
modified by the hybridization, as we have found at the end of the preceding section. With these coupling constants, 
Eq.(34) and (35) take the simple form : 



A 1 = \nTYp T -k7rTY^ 

71 ^ & 71. ^' 



(36) 



77, ^ & TJ, 1' 



(37) 



From Eq.(32) and (33) we have 



Ai _ 7l (A 2 -Ai) 

M 



|w| (|w| +7i +72) 



(38) 



A 2 
n 



A2 _7 2 (Ai-A 2 ) 
M 



|w| (|w| +71 +72) 



(39) 



where we have set 71 = T N t / N\ and 72 = T jVj / N2 ■ As expected the inverse lifetime T due to scattering within 
planes and chains has disappeared from the equations, and only are left 71 and 72 which describe physically the 
scattering from plane to chain due to hybridization. When this is carried into Eq.(36) and (37) and the summations 
on Matsubara frequency is carried out, one finds : 



k' 



A! = X[A 1X + (A 2 - Ax)— -K(p)\ - k[A 2 x + (Ai - A 2 )——K(p)} 



(40) 



A 2 = \'[A 2 x + (Ai - A 2 )J^K{p)] - k>[A lX + (A 2 - A 1 ) ] -^K(p)} (41) 

where x = In ( 1.13 u>d / T c ) and K(p) = i/>(l/2 + p ) - ip(l/2) , with p = +72)/ (27rT c ), and we have used the fact 
that in our case k / k' = 71 / 72 . Eq.(40) and (41) are essentially what one obtains for the effect of the impurities in 
a two band model [p6| . This is easy to understand physically since, because of hybridization, the in-plane scattering 
potential allows effectively an electron to scatter from plane to chain for example. However our essential point is 
that the corresponding effective inverse scattering time 7 = 71 + 72 is strongly reduced, compared to plane-plane 
or chain-chain scattering, because the hybridization is small. We note also that, in the general case, our equations 
Eq.(31) do not reduce to a simple two-band model, although we expect the physics to be similar. 

Let us briefly review the existing literature |26| for the consequences of Eq.(40) and (41) for the critical temperature, 
and add some results for the case which is of interest for us, namely the general situation with repulsive interactions. 
From Eq.(40) and (41) one gets an equation for x which is more conveniently rewritten as : 

= (x-xo){x- Xl ) 

X - X 2 

Here xo and x\ are the solutions of Eq.(9). More precisely Eq.(9) has two physically relevant solutions when AA' 
> kk'. The lower one xq gives the physical critical temperature. The higher one x\ corresponds to an unstable 
superconducting state where the order parameter has the same sign in the two bands, and therefore the interband 
coupling decreases the critical temperature with respect to the uncoupled situation instead of increasing it ( in the 
uncoupled limit xq = 1/A and xi = l/X' ; we assume for example A > A' pl[). When AA' < kk', the interband 
coupling is too strong for the state with the ill-chosen signs of the order parameter to exist and the unstable solution 
disappears ( there is still a mathematical solution with x\ < but it must be rejected physically because it does not 
satisfy the weak coupling condition x >> 1 anymore ). Finally we have ( AA' - kk' ) x 2 = ( A k + A' k' + 2 kk' )/ 
(k+k') in Eq.(42). It is easy to see that one has always xo < x 2 < x\ for AA' > kk', while x\ < x 2 < xq is satisfied 
for AA' < kk' . Now for fixed 7 , K(p) is a positive, increasing function of x + In (j/ojd) , with K(p) — > for x — > 







- oo and K(p) sa x + In (p/ /ljd) for x — > oo . Then it is easy to see graphically that, starting from xq for 7 = 0, x 
as given by Eq.(42) is an increasing function of 7 . Therefore, in all possible cases, the critical temperature decreases 
when the impurity concentration increases. 

For AA' > kk' the increase of x saturates at X2 leading to a large concentration limit T c °° for the critical temperature 
given, as found by Kusakabe p3|, by : 



1.13u D 2kk' + \k + \'k' 
T c °° ' ~ (AA' - kk')(k + k') 



(43) 



This saturation corresponds physically to the situation where strong impurity scattering has made the order parameter 
completely isotropic Ai = A2 while ( Ai - A2 ) K(p) takes an independent nonzero value. The result Eq.(43) is then 
easily rederived directly from Eq.(40) and (41). Therefore in this case the superconducting state survives very strong 
impurity scattering by becoming isotropic. This situation is analogous to the fate of very dirty anisotropic s-wave 
superconductors [El and naturally quite different from the standard behaviour of a d-wave superconductor. We note 
incidentally that, when one goes from the pure system to the very dirty one, one goes from Ai A2 < to Ai A2 > 
0, which implies that, if A + k' > A' + k for example, the chains become gapless A2 = at some specific impurity 
concentration ( given from Eq.(40) and (41) by K(p) = (k + k')/(XX' - kk') and x = (A' + k)/(XX' - kk') ). However 
this case AA' > kk' requires not too small coupling constants A and A' in both bands. In YBCO it is not clear at all 
that the CuO chains alone have a superconducting tendency. If we have A' w 0, the condition AA' > kk' cannot be 
satisfied. 

For AA' < kk' there is no saturation in the increase of x and the critical temperature goes to zero when the impurity 
concentration is increased, as found by Kusakabe |26| . However when we substitute in Eq.(42) the large x behaviour 
of K(p), we obtain : 



,_, -i n T c° _ (so - S2)M1.137/??) 
T c x 2 -xi-en(l.lfry/T°) 



(44) 



where T® is the critical temperature without impurities. Therefore T c goes to zero for a critical impurity concentration 
and we find in this case a behaviour similar to the standard d-wave result. The corresponding critical value 7 C for 7 
is given by : 



1 L13 > 

In T0 =X 2 -Xi = 



1 



2(kk' - AA') 



{X-X 1 ) +Akk' - 



(A - \')(k - k') + Akk' 



k + k' 



(45) 



The right hand side of Eq.(45) is positive and it goes from zero ( for A - A' = k - k' , or k' = ) to infinity ( for A A' = 
k k' ). Therefore the lowest possible value of 7 C is the standard d-wave result -f c = .88 T c ° , but it can easily be much 
higher since it is given by the exponential of the r.h.s. of Eq.(45). From this, one would conclude that, for a given 
7, our model is much less sensitive to impurities than standard d-waves. However when one looks at the behaviour 
of T c as a function of 7, one sees that this large value of 7 C is obtained at the end of a long tail where the value of 
T c is already quite small, as it can be seen in Fig.l. Therefore in order to fully characterize the sensitivity of T c to 
impurities, we have also to look at the initial slope found for low impurity concentration. 



In this regime we have from Eq.(42) 



_ Tc_ _ k 7 £q 
TO ~ 4 Tg x 



X2 



Xi 



4 T? [ 2 



(A - X')(k — k') +4kk' 



(46) 



2(k + k')J (X - X'Y + Akk' 



which is the result of Moskalenko and Palistrant [^6[. As they indicated, ( xq - X2 ) / ( Xq - x\ ) goes from to 1, 
which implies that the decrease of T c is always less than for the pure d-wave case. The d-wave result is obtained in 
the limiting case A - A' = k - k' or k' — > 0. On the other hand for k — > ( and fixed k' ), the critical temperature does 
not depend anymore on the impurity concentration ( this result is actually valid whatever the impurity concentration 
as it can be seen from Eq.(42)) . This limiting case is not unphysical since it corresponds to a situation where the 
density of states N c in the chains would be very small. We would still have a d-wave type order parameter with 
change of sign and nodes, but the critical temperature would be completely controlled by the planes and the chains 
would simply follow. Therefore we have the surprising result of a superconducting phase with a d-wave type order 
parameter which has a critical temperature insensitive to impurities just as an s-wave superconductor. However it 
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is fair to say that we do not expect such a limiting situation to occur in YBCO. Hence we believe that in our case 
Eq.(46) leads to a sensitivity to impurity which is always somewhat reduced with respect to the standard d-wave case, 
but not by a large factor. An example of the general behaviour of T c / Tj? as a function of our scattering rate 7 is 
given in Fig.l where we have chosen the parameters A = 1, A' = , k' = 0.5 and k takes the values 0.01 , 0.1 , 0.2 , 
0.5 , 1.5 and 5 ( k = 1.5 gives exactly an Abrikosov Gorkov law ). 

In order to conclude this section, let us now come back to the various terms we have omitted from the beginning. 
We did not take into account direct plane-chain scattering by impurities. As we mentionned, we expect physically 
this scattering to be reduced by a factor N t / N p , c compared to plane-plane or chain-chain scattering which we have 
treated. Including this plane-chain term is just what is done in the standard two-band model when the effect of 
impurities on T c is calculated. Since we have shown that our calculation lead us to the two-band model, including 
direct plane-chain scattering would just give us an additional term of the same order as the one we have found ( 
it would appear as an additional term in £0 ). Therefore our conclusions with respect to the sensitivity of T c to 
impurities are unchanged. We have also omitted mixed plane-chain terms A 12 for the order parameter. This is 
justified by our finding that, when the hybridized bands are well separated, the dominant contributions of g come 
from the block diagonal parts gn and (722 ■ Hence we have for example G\i = an gxi 021 + Qi2 .922 ^22 ■ However the 
products an &2i and a\2 a.22 are significantly different from zero only in the anticrossing region. When we calculate 
the off-diagonal self-energy, we have to sum G12 over k. Accordingly the result A12 will be smaller by a factor of 
order Nt / N PtC compared to Ai and A2, which justifies our approximation. Finally we note that going beyond the 
Born approximation will not only renormalize the impurity cross section. It will also introduce plane-chain terms £12 
in the impurity self-energy because, after scattering on a plane impurity, an electron can go to the chains through 
G12 and scatter on a chain impurity. But we have just seen that, after summation of G12 over k , one obtains a result 
which is smaller by a factor N t / N pc . This leads to a £i 2 which is smaller than £1 or £ 2 by the same factor and 
hence negligible. 



IV. DISCUSSION 



The conclusion of our study is that, in our model, the critical temperature is much less sensitive to impurities 
than it is in standard d-wave models. One basic reason is that we expect physically impurities to produce essentially 
plane-plane and chain-chain scattering. Just as in s-wave superconductors, this scattering does not affect the critical 
temperature. On the other hand T c is reduced by scattering between parts of the Fermi surface which have opposite 
signs for the order parameter, just as in standard d-wave. In our model this is due to plane-chain scattering. We have 
found that this scattering, whatever its origin, will be smaller by a factor of order t / Ef ( that is hybridization coupling 
over Fermi energy ) compared to plane-plane and chain-chain scattering. From the band structure calculations we 
expect this factor to be typically somewhere between 0.1 and 0.3 . Therefore the sensitivity of T c to impurities in our 
model is reduced by a similar factor, compared to the d-wave situation. In the specific case which we have studied 
in details and which reduces to the two-band model, we have found a further reduction of the sensitivity of T c to 
impurities with a behaviour which can vary continuously from s-wave like to d-wave like depending on the parameters. 
From our discussion of its physical origin, we expect a similar behaviour and reduction to occur in the general case. 

It is unfortunately not possible to make quantitative comparison with experiments. Indeed it seems surprisingly 
quite difficult |l^] to avoid a magnetic character for substitutional impurities in YBCO. Naturally this occurs when 
the isolated impurity atom itself has a magnetic moment. But this happens also when a non magnetic impu rity 
acquires a magnetic moment due to its interaction with the environment, as it is the case for Zn for example |p.8| . 
Naturally an impurity with a magnetic character will produce pair-breaking leading to an Abrikosov- Gorkov like law 
in any model, in agreement with what is observed experimentally. Hence these kind of experiments can not be used 
directly to eliminate theoretical models. Nevertheless, after these words of caution with respect to a simple-minded 
interpretation of impurities experimental results, we note that the reduction of T c down to 13 K for 8% Zn is obtained 
by an increase of the residual resistivity jl8) by a factor of order 10. The corresponding fi / r should be then of 
order 20 T c ( see below ). This is much more than what is necessary to destroy superconductivity within a d-wave 
model according to the Abrikosov and Gorkov law. Moreover it is useful to plot the results of Ref.17 and Ref.18 
for the critical temperature as a function of the residual resistivity. This is done on Fig. 2. It can be seen that the 
behaviour of some of our results found in Fig. 1 is quite similar to the experimental results ( we assume naturally that 
the residual resistivity is proportional to the scattering rate due to impurities ). We have naturally enough adjustable 
parameters to fit them nicely. However such a fit would be rather meaningless because, in addition to the magnetic 
problem, the scatter in the data is rather important at low impurity content and the interpretation of the data for 
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larger concentration is uncertain ( localisation effects which are not taken into account in our theory might play a 
significant role). 

If we turn to irradiation experiments, whether by electrons or by light ions, the interpretation is also not an 
easy one. It is known that high T c compounds have a critical temperature much more sensitive to irradiation than 
standard superconducting materials. However, although it is likely that most of the created defects do not have 
a magnetic character, we can not eliminate, from our knowledge on substitutional impurities, the possibility that 
some are magnetic. The experimental evidence is controversial in this respect [ p^[ . Then there is some evidence 
that localisation effects might be important since a metal - insulator transition is observed in YBCO under light ion 
irradiation, and there is no intermediate normal phase between the superconducting phase and the insulating one 
Jil]] . This sensitivity to localisation is easy to understand because of the two-dimensional nature of the CuO-i planes. 
Localisation effects arc not included in our theoretical study. It is also not clear at all that the created disorder 
can be considered as homogeneous |l2| since the resistivity measurements do not always display a sharp drop at the 
critical temperature. Similarly one might expect that the defects created by irradiation are randomly distributed at 
the microscopic scale, but it is actually quite likely that the chains are more sensitive to irradiation than the planes. 
There is finally the obvious problem of having an experimental determination of the quasiparticle lifetime produced 
by disorder. Measuring the increase in resistivity appears the best way to do it although it is far from perfect 
since, for example, it measures at best transport lifetime which we expect to be somewhat larger than quasiparticle 
lifetime. 



Notwithstanding the above problems let us try to interpret the irradiation experiments with our theoretical model, 
just to see what comes out. From Drude's law, with a typical resistivity of 100 fiQ cm at T c and a plasma frequency 
of 1.1 eV jn|, we have a typical inverse lifetime Ti / t m 2 T c . Since superconductivity disappears in d-wave when 
the inverse lifetime due to impurities is of order 2 T c we would expect that an increase of resistivity of 100 /lil cm 
leads to the suppression of T c jO]. This corresponds roughly to a decrease of T c of IK per /xfi cm (the AG law is 
essentially linear ). The experimental results of Ref.15 give a linear decrease of T c with respect to resistivity with a 
slope 0.3 K / fifl cm. In Ref.14 there is an upward curvature at low resistivity with a maximal initial slope of 0.1 K/ 
yitfi cm. Since in our model a reduction by a typical factor 1/10 with respect to the d-wave result corresponds to a 
typical choice of our parameters, we see that this last experimental result agrees with our expectation. But the result 
of Ref.14 could easily be explained for example by a larger value of the hybridization energy or by a suitable choice 
of the other parameters of our two-band model. Quite generally we have enough parameters to vary in our model 
so we can get easily agreement with these various experimental results. However it must also be kept in mind that, 
as discussed above, there are other possible physical processes which we have not taken into account and which will 
add up to produce a faster decrease of the critical temperature with the resistivity. A clear example of this is found 
in Ref.14 where a more rapid decrease is found for a set of samples and attributed to extrinsic effects, while good 
samples show an average decrease of 0.03 K/ fj.fl cm. Therefore we can consider an experimental result as an upper 
bound for our theoretical result, but it may quite well be larger than what we find. Clearly it is more difficult to 
explain in a d-wave model the slow dependence of T c on resistivity found in Ref.14 than it is in our model to explain 
the somewhat stronger dependence found in Ref.15. 

In conclusion we have seen that our model is quite coherent with the present experimental evidence. Naturally, 
with respect to this problem of the effect of impurities on the critical temperature, it would be much better to have 
experiments providing stronger constraints on theoretical models, but this might prove difficult to achieve. 

We are extremely grateful to N. Bontemps, P. Monod and J. Lesueur for very useful discussions on the experimental 
situation. 



APPENDIX 

Let us vary Eq.(9) with respect to the phonon frequencies u>d and £1 . We assume for example Sr — - Sljd / = 
- 5Q I fl < 0. The variation of Eq.(9) gives : 

Sy(2y-X* -X'*) = S(k*k'*)-Sfi*(y-X'*)-StJ,'*(y-X*) (47) 

where y = x^ 1 . The isotope effect is reduced if we show that Sy < . We set D~ x = (l + /xr)(l + fi'r ) - kk'r 2 . 
We have 6 ( k*k'* ) =2k*k'* SD / D , and Sfi* = (i* 5D / D + D ( fj.fi' - kk' ) 6r > fi* SD / D ( we assume fxfj,' - kk' 
> , which implies fiffi'* - k*k'* > ). Since we have y > A* and y > A'*, and 5D > it is enough to prove that : 
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2k*k'* - n*(y - A'*) - fi'*(y — A*) < (48) 

The left hand side is zero when the variable y is equal to Y = (2k* k'* + [i* A'* + //* + //*). We find that y\ 

< Y < yo, since we have : 

(H* + //*) 2 [(Y - A*)(Y - A'*) - k*k'*} = 
-k*k'*[\* - A'* + // - n'*} 2 - - k*k'*)[(\* - A'*) 2 + Ak*k'*} < (49) 

where yo = 1/ x o and y\ —1/ x\ . Since we have y > yo > Y , Eq.(48) is satisfied and the isotope effect is indeed 
reduced. There is no reduction only in the limiting case fifj,' = kk' and A* + ^* = A'* + fi'*. 



* Laboratoire associe au Centre National de la Recherche Scientifique et aux Universites Paris 6 et Paris 7. 
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FIG. 1. Relative variation of the critical temperature from Eq.42 as a function of our reduced effective relaxation rate 7 / 
T° , for the parameters A = 1, A' = , k' = 0.5 and for k taking the values 0.01 , 0.1 , 0.2 , 0.5 , 1.5 and 5 ( k = 1.5 gives 
exactly an AG law) . 

FIG. 2. Experimental results for the variation of the critical temperature of YBCO, due to Zn impurities, as a function of 
the residual resistivity. The filled circles are the results of Ref.17, the open squares, the crosses, the filled diamonds and the 
filled triangles are respectively the results for the batches 2, 3, 4 and 6 of Ref.18 . 
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